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Abstract 

The scattering of photons off photons at the one-loop level is investigated. We give a short review 
of the weak field limit, as given by the first order term in the series expansion of the Heisenberg- 
Euler Lagrangian. The dispersion relation for a photon in a radiation gas is presented. Based on 
this, a wave kinetic equation and a set of fluid equations is formulated. These equations describe 
the interaction between a partially coherent electromagnetic pulse and an intense radiation gas. 
The implications of the results are discussed. 
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I. INTRODUCTION 


Photon-photon scattering is a non-classical effect arising in quantum electrodynamics 
(QED) due to virtual electron-positron pairs in vacuum. Effectively, the interaction be¬ 
tween photons and these virtual pairs will result in what is known as elastic photon-photon 
collisions. Formulated as an effective field theory, using the Heisenberg-Euler Lagrangian 
3, this results in nonlinear corrections to Maxwell’s vacuum equations, which to lowest 
order in the hne structure constant are cubic in the electromagnetic (EM) held. These cor¬ 
rection takes the same form as in nonlinear optics, where the material properties of, e.g., 
optical hbres, gives rise to cubic nonlinear terms in Maxwell’s equations, so called Kerr 
nonlinearities. Since the effective self-interaction term in proportional to the hne structure 
constant squared, this means that the held strengths need to reach appreciable values until 


these ehects becomes important P|. Higher order corrections can easily be incorporated, 
as will be done here, by takii^ into account higher vertex and loop order diagrams. Pos¬ 
sible detection techniques jj, |^, as well as physical implications, such as the formation of 
light bullets Q, 0, S, I3|) photon-photon scattering have attracted a rather constant 
interest since hrst discussed (for a survey, both historical and current, of the research in 
this area, see Refs. P, |l^ |ll|, 1^ and references therein). The concept of self-trapping of 
photons due to vacuum nonlinearities has also been discussed in the context of the one¬ 


dimensional nonlinear Schrodinger equations 


Q. 


The non-trivial propagation of photons 


in strong background electromagnetic fields, due to effects of nonlinear electrodynamics, has 


been considered in a number of papers (e.g.. Ref. 


. Q). 


The main focus in these papers has 


been on the photon splitting and birefringence of vacuum, something which has attracted 
attention when it comes to the extreme magnetic helds outside magnetars ^ Q| . It has 
even been suggested that the nontrivial refractive index due to photon-photon scattering 
could induce a lensing effect in the neighbourhood of a magnetar jf^ . 

Here we will investigate the propagation of incoherent high frequency photons on a radia¬ 
tion background, both of arbitrary intensities. It will be shown that the nontrivial dispersion 
relation for a single photon on a given background gives rise to novel collective nonlinear 
effects. Furthermore, the general dispersion relation for a perturbation in the system will 
be derived, and applications will be considered. 
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II. BASIC RELATIONS 


The effective field theory of one-loop photon-photon scattering in constant background 


EM helds {E, B} can be described by the 
is the classical free field Lagrangian, and 

1 ds . 

Jq §3 ® 

■i(es)^(a^ — 6^) — 1 


/agrangian density where = —^ 


= 


[esf'ah coth(eas) cot(e6s) 


( 1 ) 


where is the electron mass, e = |e| is the electron charge, a = -l- 

b = [(^2 ^ ^2)1/2 _ ^]i/2^ ^ ^ F„feF“V4 = (^2 - E‘^)/2, ^ = FabF^^/A = -E B, 


pab ^ ^ohcdp^^l2_ 

Following Ref. 


'hus, ^ = (c? — &^)/2 and \^\ = a6. 


the above one-loop Lagrangian yields the dispersion relation for a 
test photon i n a background EM field of arbitrary field strength. This can be formulated 
according to 


1 + ^X{E^ + B^) + X -2{k-n) 

+ \{E'^ + ^2) - {k ■ Ef - {k ■ Bf = n^. 


( 2 ) 


where n = |fc|/a; is the refractive index of the background, uj is the photon frequency, k 
is the photon wavevector, and k = k/\k\. Here, we have used Abraham’s definition of the 

Tl 

momentum density of the background EM field [19|] 

E X B 


n 


n 


Furthermore, the “effective action charge” A is given by 

{d% + dl)^ 


X = 


( 3 ) 


( 4 ) 


-2d^^ + + dl)^ - 2{^d% + 

From Eqs. 0 and 0 we see that as the background EM helds vanish, so does the 
nonlinear effects, i.e., n —>• 1, as expected. 


III. DERIVATION OF THE GOVERNING EQUATIONS 

Suppose that a plane wave pulse travels through an approximately isotropic and homoge¬ 
neous medium with refractive index n. The relations E ■ B = 0 and \B\ = \E\n then holds. 
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We note that if the background medium is the plane wave held itself, the only physical solu¬ 
tion to Eq. © is n = 1, i.e., a plane wave cannot self-interact. If the background medium is 
a radiation gas, which we consider as an ensemble {Eg, Bg} of plane waves, the interaction 
contribution to the dispersion relation (j2I) will be nonzero. We hrst note that, since we are 
interested in the case of arbitrary intense helds, we have to take the effect on each photon 
in both the pulse and the gas into account. The self-interaction within the gas will also be 
nonzero, i.e., each photon in the gas will experience the refractive index n due to the gas 
and possible partially coherent plane wave pulses. Thus, taking the ensemble average over 
the radiation gas background in Eq. m, we find the relation jl 7l | 

2 + XgS 


= 

'‘'9 


1 ~ + \l^~ 2XgS’ + ^{XgS’)^ 


(5a) 


for the refractive index of a radiation gas, as experienced by a plane wave pulse. The 
relation is valid for a weakly anisotropic and inhomogeneous radiation gas, since we have 
neglected the contribution from the averaged radiation gas momentum density. Moreover, 
we have used {{kp ■ EgY) ~ {{kp ■ Bgf‘)lr? ~ ^/ 3 , where S = {\Eg\‘^), and kp denotes 
the unit wavevector of the pulse. In other words, the radiation gas is assumed to be close 
to equilibrium. The relation (|Haj) looks deceptively simple, but one has to keep in mind 
that Xg = Xging,^) is determined through Eq. (jl|). Thus, Eq. (IHafl constitutes an implicit 
relation for the refractive index. 

In order to determine Xg, we have to evaluate the derivatives of the Lagrangian (CD- For 
a radiation gas background, where each photon experiences the radiation gas through its 
refractive index, the invariants satisfy a = [(n^ — and 6 = 0. Neglecting terms 

proportional to a in the denominator of Eq. 0, we obtain the approximate expression 


for the effective action charge. Here F{ajEcnt) = (47r/a)a^(c?^ -f d'^)^c\b= 0 i ciiid Ecnt = 
TTig/e ~ 10^® V/cm is the Schwinger held j^. 

In Fig. [U we have plotted the dimensionless number XgEF,^^ as a function of the dimen¬ 
sionless variable a/Ecnt- However, since a is a function of n, Eq. (jKajl still has to be solved 
for numerically in order to obtain the refractive index as a function of the radiation gas 
energy density S. 
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FIG. 1: The dimensionless function {aE‘l^-^^/S'Ka?‘)F{a/Ecnt) x 10^ plotted as a function of the 
dimensionless variable a/Ecrit- 


The case of a photon in random motion in the background field of a plane wave pulse, 
represented by the slowly varying wave amplitude Ep, can be treated in a similar way. Taking 
the ensemble average of Eq. (ED with respect to a random ensemble of test photons yields 


nl = 


2 + XJ\Ejj 


1 - lXp\Ep\^ + ^l-2Xp\Ep\^ + \{Xp\Ep\^y 


(6a) 


Note that this can easily be generalised to incorporate the self-interaction within the gas 
(which is non-zero in general). We only need to make the replacement \Ep\^ —>■ fftot, where 
^tot = ^ Equation (IHafl has to be supplemented by an expression for the effective 

action charge, obtained from Eq. 0 according to 

OL 

Xp ~ ^^E{ap/Ecrit), (6b) 

where = [(n^ — By the replacement \Ep\'^ —>• ^tot we also include the gas 

self-interaction. 

In the weak field limit, Up^g is close to unity, and the approximation Up^g ~ 1 may therefore 
be used in the evaluation of Xg^p. We then obtain Xg^p = X = (22/45)(a/47r)E“f^. Keeping 
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the first order terms in the energy densities of the pulse and the gas, Eqs. 
yields 


Hg - 1 + lAfO*, 


and 


rip — 1 + |A|Epp, 

respectively. These results coincide with the ones obtained in Ref. j^. 


al) and (Eil 


(7a) 

(7b) 


Collective interactions 


The relation between the energy density S’ and the refractive index n{S), as given by 
Eqs. or in the weak held limit by Eq. (I7a|l . can be used to derive a wave kinetic equation, 
determining the collective dynamics of a partially coherent pulse of high frequency photons 

Q 


dl 

m 


ng{S) 


k-VI + 


|fc| dngjS) ^ = n 

ng{Sy dS ’ dk 


(8a) 


obtained using Hamilton’s ray equations for the individual photons. Here, the specihc 
intensity I{t, r, k) is normalised such that 


(|Ep|=> = iSk 


(8b) 


is the energy density of the photons, and the angular brackets is the ensemble average over 
the partially coherent photons in the pulse. 

Following Ref. the dynamics of the low frequency radiation gas, giving rise to the 
energy density S, can be formulated in terms of a coupled set of huid equations. Assuming 
a particle frame, and an equation of state P = (f/3, where P is the radiation gas pressure, 
the huid hierarchy is closed at the second moment to yield 

n \ _ S dnp(^ot) 


— V 

dt ^ V’^p(^tot 
where S{r,t) is weakly varying, and 

dn 1 


, vs = 

dt 3 UpiStot) dStot 


np{Stot) dStot dt 

S dnp(Stot) 


-VS, 


tot ? 


(9a) 


(9b) 


where ITTI is small compared to S, and Sto, = The equations (jHI) together with 

m gives the collective interaction between partially coherent high frequency photons and a 
radiation gas. 
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IV. STABILITY ANALYSIS 


We now turn our attention to the stability of perturbations of the system (jH)) and Q- 
Starting with 


/ — Jo(fc) +/i(t, r, fc), S’— S'o + S'i{t,r), 


(10a) 


and 


n = ni(t,r), 


(lOb) 


where Jq 3> |/i|, Sq S> \Si\ Sq 3> |i7i|, we derive the general dispersion relation for the 
linear perturbation. Inserting the form (HDD into Eqs. (jH)) and we obtain 


dh 


kl 


din 


^ —fc ■ VJi + ■ ^ = 0, 


dt 


UgO 


n 


ffO 


dk 


dSi 


+ -hv.(iT0 


n: 


'pO 


^0 / 
Upo^^^dt 


{S, + {\E,\^)), 


(lla) 

(llb) 


and 

(11c) 

where Ugo = ng{So), Upo = ng{Stot,o), = dng{So)/dSo, = dng{Stotfl)/dStotfl, and 


^0,11 




Jo,i d^k. 


( 12 ) 


Assuming a harmonic perturbation spectrum oc exp[z(l^ ■ r — fit)], Eqs. (Ill|l and m yields 
the general dispersion relation 


1 


-xn 




X 


1^2(1-3x)-n2ofl2(i + ^) 

ngo^ik ■K) + K‘^- 2{k ■ 
(rigofl — k ■ KY 


Jo d?k, 


(13) 


for a perturbation of the system (0 and Q. Here x = S^Sg^lnp^. 

In the case of a one-dimensional beam in the 2 ;-direction, Jq = (|Eop)5^(A;^)5(fc2 — /cq), 
and Eq. m becomes 


[liKl + JJ2)(1 - 3x) - + x)] (n,ofl - 

= ^n'g,{\Eo\^){Kl + Kj-nl,nY 

X [Kl + KYngo^-KY]. (14) 
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We can see that the characteristics between longitndinal = 0) and transverse = 0) 
pertnrbations are very different. The former is a stable perturbation, while the latter gives 
an instability growth rate 7 = —iQ according to 


2 ^ + xn^o(l^oP)] 

2’^po’^go(l + 


^ ^ 4n^o^^on^o(l + x)x(|Eop) 

I . + ^n'go{\Eo\^)Y_ 

Furthermore, if the beam has a spread in fc-space, as in the case of a Gaussian or 
Lorentzian distribution, there will be damping due to the poles of the integrand in Eq. 
m when the dimension of the physical system exceeds one. 



V. DISCUSSION AND CONCLUSIONS 


Currently, laser intensities can reach 10^^ W/cm^ and in the near future will go beyond 
this value, possibly reaching 10^^ W/cm^ j^. The intensity limit of ordinary laser tech¬ 
niques are well below the Schwinger critical field strength 2l|, but laser-plasma systems hold 
the pro mise of surpassing these limits, coming closer to the Schwinger intensity 10^® W/cm^ 


u 


at which the vacuum becomes fully nonlinear. When laser intensities reach the 


above field strengths, the plasma particles will achieve highly relativistic quiver velocities. 
Thus, th e relativ istic ponderomotive force will reach appreciable values, expelling the plasma 
particles 


n 


2 ^ 123 . Due to this expn 


sion, plasma channels will form and provide conditions for 


elastic photon-photon scattering 
will provide conditions at which Eqs 


27j |. Thus, the next generation laser-plasma systems 
j) and 0 will be applicable. 


In this paper, we have presented the equations describing the collective interaction be¬ 
tween a pulse of incoherent photons and a radiation gas close to equilibrium, at arbitrary 
intensities. The derivation of these equations was based on the one-loop Lagrangian of 
Schwinger, and the resulting refractive index of the background electromagnetic fields. The 
equations where linearised around a constant background, and the general dispersion rela¬ 
tion for the perturbation spectrum was derived. It was shown that transverse perturbations 
are unstable, and the corresponding instability growth rate was found. Applications of the 
results to the next generation laser-plasma systems were discussed. 
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